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Abstract 

We derive direct representations of the scaling functions of the 3d 0(4) model which 
are relevant for comparisons to other models, in particular QCD. This is done in 
terms of expansions in the scaling variable z = i/h^/'^. The expansions around 
z = and the corresponding asymptotic ones for z — ?■ ±oo overlap such that 
no interpolation is needed. The expansion coefficients are determined numerically 
from the data of a previous high statistics simulation of the 0(4) model on a three- 
dimensional lattice of linear extension L = 120. From the scaling function of the 
magnetization we calculate the leading asymptotic coefficients of the scaling function 
of the free energy density. As a result we obtain the universal amplitude ratio 
= 1.84(4) for the specific heat. Comparing the scaling function of the energy 
density to the data we find the non-singular part of the energy density e„s(T) with 
high precision and at the same time excellent scaling properties. 
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1 Introduction 



The aim of the paper is to provide representations of the scahng functions of the 
three-dimensional 0(4) model which can be used in tests of other models on their 
membership of the corresponding universality class. This is especially of importance 
for quantum chromodynamics (QCD) with two degenerate light-quark flavours at 
finite temperature, because it is believed [l]-[6] to belong to the 0(4) universality 
class at its chiral transition in the continuum limit. There exist already many 
parametrizations [7]- [12] of the magnetic equation of state. They differ essentially 
in the following aspects: 

a) The form of the magnetic equation of state which is initially parametrized. 
The most used form is the Widom- Griffiths (WG) form [13, 14], where both 
the scaling function and the scaling variable depend on the magnetization M . 
A second form has the advantage that the scaling variable is independent of 
M . It is therefore more appropriate for the comparison to M -data with errors. 
The two forms are completely equivalent in describing the critical behaviour 
of the model and they can be derived from each other. In principle, it is then 
only necessary to parametrize one form. 

b) The type of parametrization which is used. The parametrization has to de- 
scribe the correct general scaling laws as deduced from renormalization group 
(RG) theory, it must satisfy Griffiths's analyticity conditions [14] and take into 
account the Goldstone singularities in the low temperature phase. 

c) The input from which the parameters are determined. There are two main 
sources of information: from field-theory methods, for example e-expansions, 
small-field expansions or high temperature series etc. , and secondly, from 
Monte Carlo (MC) data. 

Our paper is inspired by the pioneering work of Toussaint [7]. In his paper the 
scaling function (of the second kind) for the order parameter was calculated based 
only on MC data and moreover these data were simulated at finite external field H . 
The main purpose of the paper was, like ours, to provide the scaling function for 
the QCD analysis. The parametrization was carried out in a third, unusual form, 
which has not been used since then, yet the Goldstone effects had still not been 
taken into account. Like the WG form, this third form has the drawback, that the 
calculation of the second form of the scaling function and in particular its derivatives 
with respect to its scaling variable is an indirect one and it is therefore unhandy. 
Moreover, scaling functions of the second kind determine via the location of their 
extrema the important pseudocritical lines. In Refs. [8] and [10] the WG form 
was used and parametrized with a combination of a low and a high temperature 
ansatz in accord with the requirements of b). The two parts were subsequently 
connected by an interpolation. The parameters were deduced exclusively from MC 
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data with finite external fields. In Refs. [9, 11] and [12] the WG form was used. 
All three papers discuss and use variants of the classical parametric representation 
of the equation of state introduced by Schofield and Josephson [15]-[17] in 1969, 
which is valid in the whole critical region. Refs. [9] and [11] differ in details of the 
parametrization and in the input: whereas [9] relies essentially on field theory input 
and uses the data for testing, Ref. [11] determines its parameters directly from fits 
to the data. In Ref. [12] the functional RG method is used to calculate the scahng 
functions. The classical representation had been invented for the WG form, at a 
time where only few MC data on small lattices were available. Today it is still used 
for all kind of calculations of universal quantities from field theory. However, as we 
shall demonstrate, it is not necessary to work with this representation. Instead we 
parametrize directly the scaling functions of the second kind, which is the preferred 
choice of the QCD community in their tests on the universality class of the 0(4) 
model. In addition, this allows us not only to make use of the magnetization data 
in the determination of the parameters, but also of the data for the susceptibilities. 

In order to broaden the tests on the universality class we calculate as well the 
scaling functions connected to the energy density and the specific heat. Here how- 
ever, the initially unknown regular or non-singular part of the energy density and/or 
specific heat is contained in the data and has to be subtracted correctly. We show 
how this can be done in principle, when the critical exponent a is negative. The 
inadequate estimate of the regular parts in a former test for 2-flavour QCD [18] 
prevented a successful outcome of the test from the beginning and led to unjustified 
conclusions. Yet, there is a derivative of the energy density, the thermal suscepti- 
bility or covariance between the energy density and the magnetization, which does 
not require a subtraction of the non-singular term. We shall use also this quantity 
for our parametrization. 

The specific model which we study here is the standard 0(4)-invariant nonhnear 
(j-model, which is defined by 

p% ^ -J 4-4 - H- ^4, (1) 

<x,y> X 

where x and y are nearest-neighbour sites on a three-dimensional hypercubic lattice, 
and c^f is a four- component unit vector at site x. The coupling J and the exter- 
nal magnetic field H are reduced quantities, that is they contain already a factor 
/3 — 1/T. In fact, we consider in the following the couphng directly as the inverse 
temperature, J = 1/T. 

It is useful to decompose the spin vector (pg into longitudinal (parallel to the 
magnetic field H) and transverse components 

4 = 4^H + 0i , With en = H/H , (2) 
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where H is the magnitude of the magnetic field. We define the energy of a spin 
configuration as 

E ^ $s4y- (3) 

<S,y> 

The lattice average 0ll of the longitudinal spin components is 

<^" = ^E4, (4) 

X 

where V = Li^ and L is the number of lattice points per direction. The partition 
function is then 

Z{T,H) = ll[d'cl>s6{$i-l)eM-^E + HVcl>\\) . (5) 

We introduce the (reduced) free energy density as usual by 

f{T,H) = -^InZ , (6) 

from which one obtains the order parameter of the system, the magnetization M, 
as 

M = = (7) 

The longitudinal susceptibility is the second derivative of — / with respect to the 
field 

XL = ^ = y((0ll^)-M^). (8) 



The energy density is 
and the specific heat 



dH 



e = ^ = (9) 

'a/3 V ' ^ ^ 



C - ^ - y{{E')-{Ef). (10) 
Finally we define the thermal susceptibility Xt as the mixed second derivative of / 

X* = ^ = (11) 

The rest of the paper is organized as follows. First we discuss the critical behaviour 
of the observables and the universal scaling functions, which we want to calculate. 

In Section 3 we describe the expansions with which we parametrize the scaling 
functions. Some details of the used simulations and the parametrizations resulting 
from the data are presented in Section 4. Here we also investigate the role of the 
non-singular terms for the scahng of the data. We close with a summary and the 
conclusions. 
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2 Critical behaviour and scaling functions 



In the thermodynamic hmit {V — )■ oo) the above defined observables show power 
law behaviour close to Tc. It is described by critical amplitudes and exponents of 
the reduced temperature t = (T — Tc)/Tc for H = and the magnetic field H for 
t = , respectively. According to RG theory the non-analytic or singular part 
of the free energy density is responsible for critical behaviour. Besides fs , the free 
energy density contains a regular or non-singular part fns- Correspondingly, the 
derivatives of fns contribute regular terms to the scaling laws, which apart from 
the cases of the energy density and the specific heat (for a < 0) are sub- leading. 
In the two-dimensional Ising model such an analytic contribution to the magnetic 
susceptibility was established [19, 20]. In Ref. [7], Toussaint makes a corresponding 
ansatz fns = ch2H^ + cjit + Cj2t^ + cjst^, which leads to an additional constant in 
Xl-i a term ~ if in M and an if-independent e„s(T). Since in our former scaling 
fits to M at Tc, e. g. in Refs. [10] and [21], we never discovered such a regular term 
we follow Privman et al. [22] and assume the non-singular part fns to have no field 
dependence, that is 

/(T, H) = fs{T, H) + fns{T) , and (12) 

e = es + ens{T) , C = C, + C„,(T) . (13) 
The regular parts do not disappear at T = Tc . We may expand e„s(T) in T at Tc 

ens{T) = e„,(Tc) + (T - T,) ■ C„,(Tc) + \iT - ■ C;^{n) + ... . (14) 

The scaling laws at if = are then for the magnetization (from now on /3 denotes 
a critical exponent) 

M = B{-tf for t < , (15) 
the longitudinal susceptibility 

XL = C+t"^ for t > , (16) 

and the energy density and the specfic heat both for f < and t > 

e = ens{T) + — rTct|t|-" , (17) 

a(l — a) 

C = Cns{T) + —\t\-'^ . (18) 
a 

For the thermal susceptibility we have for i < 

Xt = fiBT,{-tf-' . (19) 
On the critical line T = Tc or i = we have for if > the scaling laws 

M = B^H^'^ or H ^ D^M^ , (20) 
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XL = CH^'^-^ with C = B^/d . (21) 
The remaining observables scale as follows 

e = e„,(7;) + £;,//(i-")/^ , (22) 

- ^ = -i^-- , (23) 

Xt = , (24) 

where = a/A and A = is the so-called "gap exponent". 

Generalizations of these scaling laws to both non-zero t and iif-values may be derived 
from the RG scaling equation for fg 

fs(ui, U2, U3, . . .) = b-'^fsib^'Ui, b^'U2, b^'us, . . .) . (25) 

Here, the Uj with j = 1,2,... are the scaling fields, 6 is a positive scale factor and 
the yj are the RG eigenvalues. The class of our model has two relevant scaling fields 
ui — Ut, U2 — Uh with yt,yh > and infinitely many irrelevant ones with negative 
yj. The relevant scaling fields depend analytically on t and H and 

Ut = Ctt + 0{t^,H^) , Uh = ChH + 0{tH) . (26) 

The Ct,Ch are two model-dependent (positive) metric scale factors. Choosing b = 
Uf.y^'^ for H > one obtains from Eq. (25) the second form of scaling functions 
from 

fs(ut, Uh, «,>2) = «f ^Vs(«.«-^*/^^ 1, u,u-,'^/"^) . (27) 

Close to the critical point, for t, H small, Ut = Ctt, Uh = ChH , and the dependence 
on the irrelevant scaling fields becomes negligible, fs is a universal scaling function 
of Ut and Uh and 

/, = {chHY/y'^^!2{ctci''''HH-y^iy^) , (28) 

where is again a universal function. By comparison with the scaling laws one 
obtains 

yt = l/i^ , Vh = ^/t^c = A/i/ , or A = yh/yt , (29) 
and the hyperscaling relations 

2-a ^ dv, ^ p{5-l), du^p{l + 5). (30) 

Instead of working with two metric scale factors one usually introduces new tem- 
perature and field variables i = tTc/To and h = H/Hq which are chosen such that 
the scaling laws for the magnetization simplify to 

M{t = 0) = h}'^ and Ho ^ , (31) 

M{h = 0) = {-t Y and To = B'^'^T^ . (32) 
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The magnetic equation of state as derived from Eqs. (7) and (28) becomes then 

M = h'/'faiz) , (33) 

where fo is a universal scahng function with the argument 

z = i/h}'^ . (34) 

It fulfills the normalization conditions 

/g(0) = 1 , and fG{z)^-^J-zf . (35) 

Due to Eq. (7) , the corresponding scaling equation of the free energy density must 
then be 

/, = H,h'+'/'ff{z) , (36) 
where ff{z) is again a universal scaling function, and 



(37) 



Since the susceptibility xl is the derivative of M with respect to H we obtain from 
Eq. (33) 



with 



I (/g(^) - p^{z)^ . (39) 



For if — > at fixed t > 0, that is for z — > oo, the leading asymptotic term of is 
determined by Eq. (16) 

/x(^) ,=oo C^D,B'-'z-^ = R^z-^ , (40) 

where is a universal amplitude product. For ^ — )• oo the leading terms of fc and 
are identical, because for T > T^. and small magnetic field M (x H . The rest of 
our observables are related to ff{z) and /g(^) as follows 

es = -^//o/^('-"^/^/}(^) , (41) 
To 

2 



Xt 



- (^^) Hoh--'^f}{z) , (42) 
—h^^-^y^^fM . (43) 
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3 Expansions of the scaling functions 



In principle we have to parametrize only one scaling function, either ff{z) or fciz) , 
because they are related by the differential equation (37). We choose as usual faiz), 
because it is directly calculable from the magnetization data. Our representation of 
the scaling function is composed of three expansions: one around 2; = and two for 
z — )■ ±00 . In the following we derive relations between the expansion coefficients of 
ff{z) and /g(^) • We start with the expansions for small z 

CO 00 

ff{z) = Y,a^z\ fciz) = Y^hr^z^ . (44) 

n=0 n=0 

Prom Eq. (37) we obtain 



an , or a„ = — ■ . (45) 

a + n — 2 



The last equation connects the derivatives of the two scaling functions at 2; = 



/^(O) - ^^^^^-foi^). (46) 
and because of the first of the normalization conditions, Eq. (35) 

/g(0) = 6o = 1, and /^(O) = ao = ^ . (47) 

Next we consider the asymptotic expansion in the high temperature region, that is 
for z — )■ 00, or for t > and /i — > 0. Since M is an odd function of H for t > 
(Griffiths's condition), we must have 

00 

faiz) = ^-^-^d+z-^"^, (48) 

n=0 

The prefactor is the leading term of fa and = (see Eq. (40) and the remark 
after it). The corresponding ansatz for ff{z) is 

00 

//(^) = z'---J24z-'-^ , (49) 

n=0 

Using again Eq. (37) we are led to the relation 

^»«-2(^- with c-^-f , (50) 
however, the coefficient Cq is not specified by the last equation. 



7 



In the low temperature region, for t < and /i — > 0, that is for z — > — oo, massless 
Goldstone modes appear. They lead to a divergence of the transverse susceptibility 
Xt ~ . In addition also the longitudinal susceptibility xl is diverging on the 
coexistence curve. Here, the predicted divergence in three dimensions is [23]-[25] 

XL{T<n,H)r^H-'/^ . (51) 

This is equivalent to a dependence of the magnetization on if^/^ near the coexistence 
curve [8]. Therefore we make the following ansatz for foiz) in this region 

oc 

faiz) = i-zf -Y^d-i-z)--^/', (52) 

n=0 

where due to the second normalization condition in Eq. (35), dQ — 1 . The corre- 
sponding ansatz for ff{z) is 



oo 



n=0 



Inserting the two expansions into the differential equation (37) we find 

As in the high temperature phase the coefficient of the leading term, Cq , is not fixed, 
moreover c]" = 0, and C2 — —1 ■ In order to completely solve Eq. (37) for ff{z) we 
still have to find the coefficients Cq . Since a < for the 0(4) model we may proceed 
in the following way. First we consider the small 2;-expansions for 2; > 

00 °° L n i>z 

J^anz"- = AJ2 _^ = A^^-" / dyy-'Y.^^y- (55) 

n=3 n=3 + -/O ^^g 

= A^2-" /' dyy--^ [faiy) -l-hy- hy'] (56) 

= //(^) - ao - 012; - 022;^ . (57) 
That enables us to calculate Cq 

c+ = lim/^(^)^"-2 = A r dyy^-'[fG{y)-l-b,y-hy'] , (58) 

or, by partial integration 

POD 

/ dyy--' ifUy) - /^(O) - yfm] ■ (59) 
Jo 

In the same manner we can calculate Cg from an integral over negative z by starting 
from 

ff{z) = ao + aiz + a2z'' + A{-zf-" f' dyi-y)^-' [faiy) - 1 - hy - b2y'] , (60) 

J z 



A 

Cn = 



8 



and taking the limit 2; — > — oo 



Co 



lim ffiz){-zr-' = A r dy{-yr-' [foiy) - 1 - hy - b2y'] , (61) 

J —00 



or 



Co = ^ / dy{-yr-' [/^(y) - /^(O) - y/^(0)] . (62) 

The function ff{z) is universal as a whole and so are each of its expansion coefficients, 
that is Cq and Cq are universal. In fact, it can be shown, that these coefficients are 
the universal products of critical amplitudes 

4 = ff{B^YB-^^^'^ , (63) 

where 

ff = T , (64) 

—a[l — a)(2 — a) 

are the critical amplitudes of the free energy density for H = and t ^ 

f = fns(T) + f^\t\'-- . (65) 
Because of Eqs. (63) and (64) we can now calculate the universal ratio 

A+ ct 



A- ' c, 



(66) 







from the two integrals in Eqs. (59) and (62). A similar formula is known for the 
magnetic equation of state in the Widom-Griffiths form [26]. It was for example 
used to determine the ratio A'^ /A~ in the case of the 0(2) model [27]. 



4 The parametrizations 



The data we use in the following to parametrize the scaling functions were all ob- 
tained from simulations described in detail in Ref. [21]. We repeat here only the 
main features of these simulations. They were performed on three-dimensional lat- 
tices with periodic boundary conditions and linear extension L = 120. The coupling 
constant region which was explored is 0.90 < J < 1.2, the magnetic field was varied 
from H = 0.0001 to H = 0.007. In general 100000 measurements were done at 
each fixed H and J. We have reevaluated the raw data to obtain the magnetiza- 
tion, the longitudinal and the thermal susceptibilities, the energy density and the 
specific heat. Due to the large spatial volume of the lattice, most of the finite size 
effects have disappeared from the data. This is in particular true for the energy 
density and the magnetization, to a smaller extent also for the susceptibilities and 
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the specific heat. We shall discuss the remaining effects when the scahng of the re- 
spective observables is analyzed. A further source of difficulties is, at larger \t\ and 
if- values, the possible appearance of corrections to scaling because of the influence 
of irrelevant scaling fields. These violations of scaling should be visible in the scaling 
plots for the data. However, as we shall show below, we find perfect scaling prop- 
erties without any sign of these corrections to scahng for our values of t and H for 
the energy density and the magnetization. This not the case for the corresponding 
non-linear 0(2) model, see e. g. Fig. 5 of Ref. [28], where strong scaling violations 
were found for T < Tf.. An explanation for this striking difference can be inferred 
from two papers, by Hasenbusch and Torok [29] for = 2 and by Hasenbusch [30] 
for A" = 4. In these papers the leading corrections to scaling could be ehminated 
by using instead of the non-linear 0(A") invariant models the corresponding 0(A") 
symmetric 0^ models and the tuning of the additional parameter A . It turned out 
that the optimal parameter value is 2.1 for = 2 and 12.5(4.0) for A^ = 4. The 
non- linear case corresponds to A = oo. One expects therefore to find significantly 
weaker corrections for A^ = 4 as compared to A^ = 2. The remarkable lack of scahng 
corrections had already been noted in Ref. [8], where for the first time the critical 
exponents for the 3d 0(4) model were determined from magnetization data at finite 
external fields. Later, in Ref. [10], the result for the exponent S could be improved 
with better data and fits where correction-to-scaling terms had been taken into ac- 
count. Yet, these terms were contributing at best marginally and could as well be 
neglected. 

In order to define our variables t, i, h and z we use the same critical amplitudes, 
temperature and exponent values as in Ref. [10] and [21]. These are 

Jc = T"^ = 0.93590 , To = 1.093 , Hq = 4.845 , (67) 

p = 0.380 , 5 = 4.824 , A = 1.83312 , (68) 
and from the hyperscaling relations 

a = -0.2131 , ly = 0.7377 , 7 = 1.4531 . (69) 

We have compared the exponents from Eqs. (68) and (69) to the field theory results 
of Guida and Zinn- Justin [31], displayed in their Table 3. Apart from the value 
for S = 4.824(9) (from [10]) which corresponds to the value 77 = 0.0302(16) via 
the hyperscaling relation t] = {5 — S)/{1 + 6) our numbers are always close to the 
respective central values of Guida and Zinn- Justin and they are always inside their 
error bars. 

We procedc in the following way. First we calulate fdz) from the magnetization 
data and fit the large z parts to the asymptotic expansions. The small z region is 
more intricate, because the derivative —fc{z) has a peak for 2; > 0, which determines 
the pseudocritical line. In order to model the corresponding variation properly we 
fit directly the derivative for small z. It can be obtained either from the data for 
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Figure 1: The scaling function fc = Mh~^/^ as a function of 2; = th~^^^. The 
dashed hues show the asymptotic expansions, the sohd hne the Taylor expansion at 
2; = 0. The numbers refer to the different J = 1/T- values of the data. 



Xt, Eq. (43), or from xl and M, using Eqs. (33), (38) and (39). After completion of 
the parametrization of fciz) we compute the leading asymptotic coefficients Cg of 
ff{z). The scaling function of the free energy density is then also entirely known. 
In the next step we determine the non-singular contributions to the energy density 
and the specific heat, first at Tc and then at all our T- values. We show that the 
results for e^s shape a smooth function of T, where is not a distinguished point. 
With this function e„s(T) it is then possible to compare the scaling functions to the 
data for the energy density and specific heat. 

In Fig. 1 we show the data obtained from the magnetization for the scaling 
function fciz) and our parametrizations. Obviously the data scale very well, apart 
from the data for H = 0.0001 and z < 0, which show some finite size effect (in the 
figure at 2; = —1.534 and —5.219). We have fitted fc in the asymptotic regions 
with the first three terms of the respective expansions from Eqs. (48) and (52). For 
the positive 2;-range [1.5, 15] we found the coefficients 

4 = 1.10599 ± 0.00555 , df = -1.31829 ± 0.1087 , 4 = 1.5884 ± 0.4646 . (70) 

In the negative z-range [—10, —1] we discarded the data with H = 0.0001 and 
obtained 

4 = 1, 4 = 0.273651 ± 0.002933 , 4 = 0.0036058 ± 0.004875 . (71) 
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Figure 2: The derivative —fciz) as a function oi z = ih^^^^. The filled circles denote 
the data calculated from Xt, the crosses the data obtained from xl and M . The 
dashed lines show the asymptotic expansions, the solid line the Taylor expansion 
around z = 0. 

Since (Iq = we have a new value for this quantity, which is compatible with 
the old values = 1.084(18) from Ref. [10] and R^ = 1.12(11) from Ref. [9] but 
somewhat more accurate. Astonishingly, the asymptotic expansions describe the 
function /g(^) very well down to rather small |2;|-values, and as can be seen in 
Fig. 1 they overlap with our approximation to the Taylor expansion at 2; = . As 
mentioned already, we use the derivative —JqIz) to determine the coefficients of the 
Taylor expansion. In Fig. 2 we show the data which we obtained from Xt (filled 
circles) and xl and M (crosses) for the derivative. Obviously, the data involving 
Xl suffer from large finite size effects in the whole low temperature region {z < 0) 
for already moderately small if-values. This behaviour is known and a consequence 
of Eq. (51), the divergence of xl near the coexistence line. In contrast to that, 
the data from xt show a consistent scaling behaviour for z < (apart from the 
H = 0.0001 point for J = 0.94 at ^ = -1.534). For z > but close to the critical 
point we find still larger finite size effects for the Xi-data as compared to the ones 
for the Xt-data, because of the stronger divergence of xl on the critical line. At 
larger positive z-values beyond the peak region we observe however, that the Xt- 
data are systematically smaller than the XL-data. The reason for this lies in the 
cluster update [32] which was used to produce the data. That update diminishes 
very efficiently the autocorrelation time for the order parameter, but is less efficient 
for the energy density. For increasing temperature and/or decreasing H the cluster 
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Figure 3: The quantity —fciO) as a function of the magnetic field H. The filled 
data points were used for the fit, the solid line shows the fit result. 



size drops and if the number of cluster updates is not correspondingly increased the 
autocorrelation times for the energy density data increase faster than those for the 
magnetization data. An increasing decorrelation of e and M takes then place with 
increasing z, that is the Xf-data become too small. In view of all these considerations 
we use in our Taylor fits only xrdata for z < 1.3 and in the 2;-interval [1.3, 2] both 
types of data. 

From our data at Tc we have calculated an additional data point at z = 0, 
denoted by a star in Fig. 2. To this end we have used Eq. (43) at Tc 

-fciO) = ^,h-^'-'^/\,{n) . (72) 

c 

The corresponding data are shown in Fig. 3. At small H we have again finite size 
effects, at large H corrections to scaling set in. A fit to the remaining data points 
leads to the result — /g'(O) = 0.3173(5). We have approximated the Taylor series 
with polynomials. Instead of using a single very high order polynomial, we started 
with two fits to 6. order polynomials, that is for JqIz) with the ansatz 

f^{z) = h + 2b2Z + 3hz^ + 464^^ + 565/ + 666^^ , (73) 

in the overlapping ^-intervals [—2.5,0.75] and [—0.75,2]. The lowest coefficients 
of the two polynomials should, at the end of the process, coincide to generate a 
smooth parametrization close to 2; = 0. In the second step we took therefore the 
arithmetic averages of the results of the two fits for bi and 62, fixed them and repeated 
the two fits to determine the remaining coefficients. In step three 63 was fixed by 
averaging the corresponding means of the first and the second step. The last fits 
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were performed with fixed 61,62 and 63 in the intervals [—2.5,0.75] and [—0.5,2], 
including some points from the fits to the asymptotic regions. Our final result is 

bo = l , bi = -0.3166125 ± 0.000534 , (74) 
62 = -0.04112553 ± 0.001290 , 63 = 0.00384019 ± 0.000667 . (75) 

The remaining coefficients are different for 2; < and 2; > 0. We find for 2; > 

6| = 0.006705475 d= 0.001704 , 6+ = 0.0047342 d= 0.001429 , (76) 

b+ = -0.001931267 ± 0.000312 , (77) 

and for 2; < 

6J = 0.007100450 ± 0.000160 , 65 = 0.0023729 ± 0.000095 , (78) 

feg = 0.000272312 ± 0.000021 . (79) 

We note that 6j and still coincide inside their error bars. In Fig. 2 we have 
plotted the respective approximations to —foiz) in the 2;-ranges [—2.5, 0] and [0, 2.5]. 
Obviously, there is a large range for 2; < and a shorter range for 2; > where the 
approximations overlap and coincide with the respective asymptotic expansions. In 
the rest of the paper we use therefore the Taylor expansions in the 2;-range [—2, 1.95] 
and outside the asymptotic expressions. 

It is now straightforward to calculate the coefficients of the leading asymptotic 
terms of ff{z) from Eqs. (59) and (62). We find 

c+ = 0.422059886 ± 0.010595 , Cg = 0.229176194 ± 0.010669 . (80) 

The errors of have been determined using the complete correlation matrix of 
the contributing parameters. The main contributions to are coming from the 
two terms proportional to bi and 62. The second, larger term is the same for both 
coefficients, the first only changes the sign, that is for we have the sum, for Cq 
the difference of these terms and as a consequence the arc strongly correlated. 
That allows us to estimate the correlation between Cq and Cq to C+_ = —af + cig, 
where the ai^2 are the errors of the two terms. From Eqs. (66) and (80) we obtain 
then the universal ratio 

— = 1.842 ±0.043 . (81) 

Our value for the ratio is in agreement with the final estimates found in Refs. [9], 
1.91(10), and [11], 1-8(2) . We note however, that our error estimate does not include 
a possible variation of the critical exponents used. 

With the results from Eq. (80) we have completely specified the parametrizations 
of fdz) and of ff{z). In order to test the scaling function of the free energy density 
and/or its temperature derivatives we still need the regular contributions. Quite 
generally, the exponent a must be always less than 1 for continuous transitions. 



14 



because the energy density at Tc is finite and there is no latent heat (which would 
be possible for a = 1). From Eqs. (17) and (22) we know therefore that 



e,(T„if = 0) = , or e(T„ i7 = 0) = e„,(T,) 



f82) 



We may then determine e„s(Tc) from our data for the energy density at T^. In Fig. 4 
we see that the data at T = fulfill the expected scaling law, Eq. (22) and Eq. 
(41) accordingly 

6.(T,) = -S//„/,(i-)/Aj^(o) . (83) 

Moreover, there is no sign of a finite size dependence at small external fields. We 
have fitted the data directly with Eq. (22) and find 

e„^(Tj = -0.991888(13) , = -0.8500(06) , /}(0) = 0.47723(32) . (84) 

Evidently we have a very precise result for the non-singular part of the energy density 
at Tc. The derivative //(O) can also be calculated from bi = /^(O) and Eq. (46). 
That leads to //(O) = 0.47843(81) and is consistent with the previous result. 

For positive a the specific heat is diverging and it is therefore unclear how to 
determine the regular term Cns{Tc). If however a is negative, then Cs — "ies/T 
disappears at Tc for h = 0, because of Eq. (23), the finite regular term remains and 
we can calculate C„s(Tc). Usually, the specific heat has nevertheless a sharp peak 



-0 98 



-0 . 99 



-1 00 



-1,01 



-1 02 



-1 . 03 




, 005 



0,01 



015 



Figure 4: The energy density at (open circles) function of /^(i-")/^ . The filled 
circle shows the extrapolation to h = 0, that is e„s(T'c), the points are connected by 
straight lines. 
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at Tc and H = 0, which just means that the critical amphtudes are negative. That 
sign is taken care of by the factors 1/a and 1/ac, respectively in our amplitude 
definitions, so that and are positive. In Fig. 5 we show our data for C — 2e/T 
at Tc and the corresponding fit to Eqs. (23) and (42). We find 

(Cns - (Tc) = 6.2669(195) , and C„,(Tc) = 4.4103(195) , (85) 

Ac = 0.32041(468) , or /}'(0) = 0.7151(104) , (86) 
which is compatible to the value ff{0) = 0.7075(221) calculated from 62 ■ 

In Fig. 6 we show the scaling function ff{z) which is obtained from our parametri- 
zation. We plot the asymptotic expansions and the Taylor expansion separately and 
find that they are overlapping in the same regions as for fciz), that is our calculation 
of the is consistent. The data which we also show have been calculated assuming 
that 

ensiT) = eUTc) + (T - T,) ■ C„,(T,) , (87) 

where the numbers are from our fits at T = Tc. We note here, that without the term 
proportional to Cns(Tc) the data would scale nowhere apart from the point z = 0. 
In Fig. 6 we observe scaling for small \z\ and coincidence with the predicted scaling 
function, but not at larger values. The reason for that is the assumption, Eq. (87), 
for the function e„s(T), which leads to inaccurate e„s-values for larger |T — Tc|. We 
can test this and the scaling behaviour at fixed T and varying H at the same time. 
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10 




Figure 6: The scaling function ff{z) as a function of 2; = th~^/^. The dashed hues 
show the asymptotic expansions, the sohd hne the polynomial approximations for 
small z. The star at 2; = is the result from Eq. (84). 



If we have scaling then Eq. (41) must hold and 

eUT) = e(r, H) + T^^h^^-)/^ f {z) , (88) 

-1-0 

where the e(T, H) are the energy density data at fixed T and ff{z) is the predicted 
scaling function. The test is successful, if we obtain the same value ens{T) for all H 
inside the error bars. This is indeed the case. The errors of the averages are tiny. 
They vary for T < between 6 ■ 10~® and 1.4 ■ 10~^, and in the hot phase they 
increase to 3.5 ■ 10~^. In Fig. 7 we compare the found values for e„s(T) with the 
approximation from Eq. (87), where we have used our previously calculated numbers 
from Eqs. (84) and (85). Though the differences in Fig. 7 do not seem to be large, 
they are the reason for the deviation of the data from the scaling function ff{z) in 
Fig. 6. We have fitted our results for ens(T) with the Taylor expansion, Eq. (14), 
up to the third derivative of e„s and find 

ens{T,) = -0.991792(28) , C„,(r,) = 4.3910(14) , (89) 
C;^(T,) = 8.448(108), Cl(T,) = 42.79 ± 5.13 . (90) 

The new result for Cns{Tc) is in agreement with the one we obtained from our 
admittedly far extrapolation of C — 2e/T in h a.t Tc . We have used the results 
from the Taylor expansion fit to approximate the functions e„s(T) and C„s(T) in 
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Figure 7: The non-singular part of the energy density e„s(T) as a function of T — Tc 
(filled diamonds). The dashed line shows the approximation from Eq. (87). 

the calculation of the scaling functions and /" from our data. As can be seen 
from Fig. 8 we find now perfect scaling for the energy density, even for relatively 
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Figure 8: The scaling function f'fiz) as a function of z = th The line shows 

our parametrization, the data have been calculated using Eqs. (89) and (90). 
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Figure 9: The scaling function /j( function oi z = th The hne shows 

our parametrization, the data have been calculated using Eqs. (89) and (90), the 
star at 2; = is the result from Eq. (86). 



large |T — Tc| -values. In Fig. 9 we compare our parametrization for /j-' with the data. 
Since the specific heat is proportional to the fluctuation of the energy density, the 
data are not as precise as for /j, especially for high temperatures. Nevertheless we 
observe satisfactory scaling and a further confirmation for our parametrization. 

Finally we show in Fig. 10 the third derivative of the scaling function ff{z) with 
respect to z. It controls the singular behaviour of the third derivative of the free 
energy density with respect to temperature, 

and is the first thermal derivative of the free energy density that diverges at Tc in 
the limit if — )■ 0. It thus allows a discussion of critical behaviour resulting from 
the structure of the singular part of the free energy density without the need of 
determining non-singular contributions to the free energy density. This property of 
the three-dimensional 0(4) model has been exploited in the discussion of critical 
behaviour in the vicinity of the chiral phase transition of 2-flavor QCD [33, 34]. 
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5 Summary and conclusions 



In this paper we have investigated the scahng functions of the three-dimensional 
0(4) model, which can be derived from the singular part of the free energy density. 
In contrast to other papers [8]- [12] where the scaling functions are parametrized 
in the Widom-Griffiths form, we have chosen a form which is preferred in tests 
of other models such as QCD. Here the scaling variable z is independent of the 
magnetization. These scaling functions can of course be derived from those of the 
other type, however the explicit functional dependence on z is of value. The major 
advantages of our parametrization are 

1) derivatives with respect to the variable z can be taken directly and not im- 
plicitly, the positions of extrema are easily calculable. 

2) the parameters are determined from direct fits to the data, one can immedi- 
ately judge how well the data are represented. 

In order to carry out the parametrization we have used the best presently available 
data set for finite magnetic fields [21]. Furthermore, we have tackled the problem 
of the scaling of the energy density and the specific heat in the 3(i 0(4) model, to 
our knowledge for the first time at all, and we were able to clarify the role of the 
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non-singular part of the energy density for scaling. This is very important for all 
corresponding checks of QCD with two light flavours (see for example Ref. [18]). 
Our approach in some more detail was the following: 

We have parametrized the scaling function fdz) of the magnetization with asymp- 
totic expansions for z — > ±00 and Taylor expansions around z — 0. The knowl- 
edge of the expansion coefficients of fdz) enabled us to derive the corresponding 
coefficients for the scaling function ff{z) of the singular part of the free energy 
density. In particular, we could calculate the leading asymptotic coeffcients Cg of 
ff{z) and thereby determine the universal amplitude ratio for the specific heat to 
A'^/A~ — 1.842(43) . In the following we have tested our data for the energy density 
and the specific heat with the respective scaling functions. To this end we have de- 
termined the non-singular parts of the two observables at Tc . Whereas this is always 
possible for the energy density, we could do that for the specific heat only, because 
in the three-dimensional 0(4) model the exponent a is negative. With these results 
we found scahng for the energy density in the neighbourhood of Tc or small \z\ but 
not outside. As it turned out, we can achieve scaling for all our T or J- values, if 
we use the correct values of e„s(T). The latter can be calculated from our function 
ff{z) if the data show scaling in h. We found indeed perfect scaling in h and very 
accurate e„s(T)-values. They form a smooth function of T in accord with the results 
from Tc . As expected, Tc is not a distinguished point of this function. For other 
models the regular parts of the energy density are of course different and may be 
not so easy calculable, if for example correction-to-scaling terms spoil /^.-scaling. A 
test of the energy density and/or specific heat on scaling may then be problematic. 

As a last result of our parametrization we quote the peak positions of scaling 
functions which are relevant for the definition of pseudocritical lines. In all cases 
the peaks are very fiat and symmetric around the peak positions Zp. We found 
the value = 1.374(30) for in -/^ the value z^'^ = 0.74(4) and for -/}' the 
value Zp'^ — —0.38(8). The corresponding observables are xl, Xt, and C — 2e/T, the 
upper index of Zp denotes the number of derivatives of the free energy density with 
respect to T and H . Evidently, the pseudocritical lines can be rather different. If, 
for example, the peak position of Xt is used to define the pseudocritical temperature, 
it will be closer to Tc as in the case of xl- 
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